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Abstract. We investigate the effects of rare regions on the dynamics of Ising 
magnets with planar defects, i.e., disorder perfectly correlated in two dimensions. 
In these systems, the magnetic phase transition is smeared because static long- 
range order can develop on isolated rare regions. We first study an infinite-range 
model by numerically solving local dynamic mean-field equations. Then we use 
extremal statistics and scaling arguments to discuss the dynamics beyond mean- 
field theory. In the tail region of the smeared transition the dynamics is even 
slower than in a conventional Griffiths phase: the spin autocorrelation function 
decays like a stretched exponential at intermediate times before approaching the 
exponentially small equilibrium value following a power law at late times. 
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1. Introduction 

In recent years, the influence of quenched disorder on phase transitions and critical 
points has reattracted considerable attention. A particularly interesting source of new 
phenomena are the effects of rare spatial regions which can be locally in the wrong 
phase. These rare regions or "Griffiths islands" fluctuate very slowly because flipping 
them requires changing the order parameter in a large volume. Griffiths [1] showed 
that this leads to a non-analytic free energy everywhere in an entire temperature region 
close to the transition, which is now known as the Griffiths phase [2] or the Griffiths 
region. In generic classical systems with uncorrelated disorder, the contribution of 
the Griffiths singularities to thermodynamic (equilibrium) observables is very weak 
since the singularity in the free energy is only an essential one [1,3]. In contrast, the 
consequences for the dynamics are much more severe with the rare regions dominating 
the behavior for long times. In the Griffiths region, the spin autocorrelation function 
C(t) decays very slowly with time t, like lnC(t) ~ — (lnt) d ^ d ^ 1 " > for Ising systems 
[2,4-7], and like lnC(t) ~ — i 1 / 2 for Heisenberg systems [7,8]. Many of these results 
have been confirmed by rigorous calculations of disordered Ising models [9-11]. 

The effects of impurities are greatly enhanced by long-range spatial disorder 
correlations. In an Ising model with linear defects, the thermodynamic Griffiths 
singularities are given by power-laws with the average susceptibility actually diverging 
in a finite temperature region. The critical point itself is of exotic infinite-randomness 
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type and displays activated rather than power-law scaling. This was first found in 
the McCoy- Wu model [12, 13], a disordered 2D Ising model in which the disorder is 
perfectly correlated in one dimension and uncorrelated in the other. Later it was 
studied in great detail in the context of the quantum phase transition of the random 
transverse-field Ising model where the imaginary time dimensions plays the role of the 
"correlated" direction [14]. 

Recently, it has been shown that even stronger effects than these power-law 
(quantum) Griffiths singularities can occur in Ising models with planar defects 
[15,16]. Because the disorder is perfectly correlated in two directions, the effective 
dimensionality of the rare regions is two. Therefore, an isolated rare region can undergo 
the phase transition independently from the bulk system. This leads to a destruction 
of the global sharp phase transition by smearing. Similar smearing effects have also 
been found in itinerant quantum magnets [17] and at a non-equilibrium transition in 
the presence of linear defects [18]. In all of these cases, the effective dimensionality 
of a rare region is above the lower critical dimension of the problem which allows the 
rare region to order independently. 

In this paper, we investigate the dynamics of an Ising model with planar defects 
in the vicinity of this smeared phase transition. We show that the dynamics in the 
"tail" of the smeared transition is even slower than in a conventional Griffiths phase. 
The autocorrelation function decays like a stretched exponential at intermediate times 
followed by power-law behavior at late times. The paper is organized as follows. 
In section we introduce the model and briefly summarize the results for the 
thermodynamics of the smeared transition to the extent necessary for the discussion 
of the dynamics. In sectional we numerically solve dynamical mean- field equations 
and compare the results to predictions based on a variational approach. In section 0] 
we use general scaling arguments in connection with extremal statistics to go beyond 
mean- field theory. Conclusions are presented in section [5j 

2. The model and its thermodynamics 

For definiteness we consider a (/-dimensional Ising model with bond disorder 
completely correlated in the dc = 2 directions x\ and xi but uncorrelated in the 
remaining d± = d — dc directions X3, . . . ,Xd- Classical Ising spins S x = ±1 reside on 
the sites x of a hypercubic lattice. They interact via nearest-neighbor interactions. In 
the clean system all interactions are independent of the lattice site; their values are Jo 
for bonds in the uncorrelated directions and Jo/4 for bonds in the x\ and X2 directions. 
We have chosen this scaling of the interactions because it leads to mean-field equations 
compatible with the corresponding thermodynamic mean-field theory [15]. The planar 
defects are modelled by weakening the bonds in (x\, a^-planes at random positions 
in the uncorrelated directions. The system effectively consists of blocks separated by 
parallel planes of weak bonds. The Hamiltonian of the system is given by: 



where &j is a unit vector in direction Xj, x_l is the projection of x on the uncorrelated 
directions £3, . . . , Xd, and J(xj_) is the random coupling constant in the (x± , X2)-planes. 
The values of J(xj_) are drawn from a binary distribution: 




(1) 



X 



x j=3 



J(x±) = 



c Jo with probability p 
Jo with probability (1 — p) 



(2) 
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characterized by the concentration p and the relative strength c of the weak bonds 
(0 < p < 1,0 < c < 1). 

We now briefly summarize the results on the thermodynamics of this system 
[15, 16] to the extent necessary for the discussion of the dynamics. The clean system 
(p = 0) undergoes a magnetic phase transition at some temperature T®, with the order 
parameter being the total magnetization 

m= ^E^), (3) 

X 

where V is the system volume, and (• • •} is the thermodynamic average. In the 
disordered system (p > 0), a crucial role is played by rare strong disorder fluctuations: 
There is a small but non-zero probability for finding large spatial regions in xj_ 
direction devoid of any impurities. These regions can be locally in the ordered phase 
even if the bulk system is in the disordered phase. Because of the disorder correlations, 
these rare regions are infinite in the x\ and xi directions but finite in the remaining 
directions. Thus, each rare region is equivalent to a two-dimensional Ising system 
and can undergo a real phase transition independently of the rest of the system. The 
resulting effect is much stronger than the conventional Griffiths effects [1,2]: The global 
phase transition is destroyed by smearing [15]; and the order parameter develops very 
inhomogeneously in space with different parts of the system (i.e., different xj_ regions) 
ordering independently at different temperatures. Correspondingly, the correlation 
length in x^ direction remains finite for all temperatures. 

Using extremal statistics, the leading thermodynamic behavior can be easily 
worked out in the 'tail' of the smeared transition, i.e., in the parameter region where 
a few islands have developed static order but their density is so small that they can 
be treated as independent. The probability w of finding a large region of linear size 
Lr (in x^-space) devoid of any impurities is, up to pre-exponential factors, given by 

w ~ (1 - p) L ^ = exp[L^ ln(l - p)} . (4) 

As discussed above, such a rare region develops static long-range (ferromagnetic) order 
at some temperature T c (Lr) below the clean bulk critical temperature T c °. The value 
of T c (Ln) varies with the size of the rare region; the largest islands will develop long- 
rage order closest to the clean critical point. According to finite size scaling we obtain 

T c {L)-T° c =r c {L R ) = -AL-\ (5) 

where <fi is the finite-size scaling shift exponent of the clean system and A is the 
amplitude for the crossover from d dimensions to a slab geometry infinite in two 
dimensions but finite in d± — d — 2 dimensions. The reduced temperature r = T — T c ° 
measures the distance form the clean critical point. If hyperscaling is valid, the finite- 
size shift exponent fulfills (f> = \ jv. Combining Q and JSJ we obtain the probability 
for finding a rare region that becomes critical at r c as 

w{r c ) ~ exp(-B \r c \- d±/,p ) (for r c -»• 0-) (6) 

where the constant B is given by B — — log(l — p)A d± ^. The total (or average) 
order parameter m is obtained by integrating over all rare regions which are ordered 
at r, i.e., all rare regions having r c > r. Since the functional dependence on r of 
the order parameter on a given island is of power-law type it does not enter the 
leading exponentials but only the pre-exponential factors. To exponential accuracy, 
we therefore obtain 

m - exp(-B \r\- d±/ ' 1 ') (for r -»■ 0-) . (7) 
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Thus, the total magnetization develops an exponential tail towards the disordered 
phase which reaches all the way to clean critical point. Analogous estimates show 
that the homogeneous magnetic susceptibility does not diverge anywhere in the tail 
region of the smeared transition because the exponentially decreasing island density 
overcomes the power-law divergence of the susceptibility of an individual island. 
However, there is an essential singularity at the clean critical point produced by the 
vanishing density of ordered islands. 

3. Local dynamic mean-field approach 

We now turn to the main topic of this paper, viz., the dynamic behavior in the vicinity 
of the smeared phase transition. We consider a purely relaxational local dynamics 
without any conservation laws, i.e., model A in the Hohenberg-Halperin classification 
[19]. Microscopically, it can be realized, e.g., by the Glauber algorithm [20]. 

In this section, we study the dynamics of the Ising model Q by numerically 
solving local dynamic mean-field equations. For definiteness, we consider d = 3 spatial 
dimensions, i.e., d± = 1 and dc = 2. A local dynamic mean-field theory of the 
Hamiltonian can be derived, e.g., by following Ref. [21]. The plane magnetizations 



where Lq is the system size in correlated direction, then fulfill the local mean-field 
equations of motion 



Here, a is a microscopic relaxation rate which determines the overall time scale. 
Alternatively, the same equations can be obtained without further approximation from 
the infinite-range Hamiltonian used in the thermodynamic mean- field theory [15]. 

Note that the local mean-field equations contain the full spatial dependence of 
the magnetization on the uncorrelated dimension. In contrast to a global mean- 
field theory these local equations are therefore capable of describing rare region 
physics. The stationary solution of J3J corresponds to the thermodynamic local mean- 
field theory studied in Ref. [15]. The clean model (p = 0) undergoes a magnetic 
phase transition at Xj? = 3 Jo, displaying mean-field critical behavior. Below Tj?, 
the diluted system (p > 0) develops an exponential magnetization tail of the form 
m ~ exp(- J B|T c ° - T^ 1 / 2 ) following the prediction 10 with d± = 1 and = 2. 

To study the dynamics, we have numerically integrated the equations Q for 
a = 1, Jo = 1, c = 0.2, dilutions of p = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7and 0.8 and system sizes 
(uncorrelated direction) up to L± = 10 6 averaging over up to 25 disorder realizations. 
Let us first analyze the behavior at the clean critical temperature TJ? = 3.0 (at the 
threshold of rare region ordering). The left panel of figure ^ shows the time evolution 
of the total magnetization m = (1/L±) m x at Tj? = 3 for several dilutions. The 
figure shows that the long-time decay of the magnetization is governed by a stretched 
exponential of the form lnm(i) = —Et^ 1 ^. The decay constant E increases with 
increasing dilution. The right panel of figure ^ shows £ as a function of the dilution 
p. Clearly, the decay constant depends linearly on [— ln(l — p)] 2 ' 3 . 

This behavior can be understood using variational arguments similar to that of 
Bray [6,7] for conventional Griffiths effects. At the clean critical point, the correlation 




(8) 



— -r, m 'x = -m x + tanh [(Joto^-i 



a dt 



+ J{x)m x + J m x+ i)/T] 



(9) 
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Figure 1. Left: Time evolution of the magnetization at the clean critical 
temperature T® = 3 for several impurity concentrations p starting from a fully 
polarized state at to = to t max = 2000. Parameters are Lj_ = 10 5 , Jo = 1, 
c = 0.2. The data represent averages of 25 disorder configurations. The long-time 
behavior of m follows In m(t) ~ —Et^ 1 / 3 . Error bars representing the standard 
deviation of m are shown for p = 0.3. Right: Decay constant E of the stretched 
exponential decay of the magnetization as a function of the impurity concentration 
p. The solid line is a linear fit to [-ln(l -p)] 2 / 3 . 



time £t of an impurity- free island depends on its linear size L R via a power law £t ~ L Z R . 
Here, z is the clean dynamical exponent whose value is 2 for our infinite-range model. 
The rare region contribution to the time-dependent magnetization m(t) can obtained 
by summing over the exponential time dependencies of the individual islands with 
Li?-dependent correlation time. Using Q we obtain to leading exponential accuracy 

m(i) ~ J dL R exp [L R bx(l - p) - t/£ t (L R )] . (10) 

This integral can easily be evaluated variationally, i.e., within the saddle-point 
method. The main contribution to m at time t comes from islands of size L R p (t) ~ 
[t/(— ln(l — p))] 1 / 3 . The leading long-time decay of the magnetization at the clean 
critical point is then given by a stretched exponential 

lnm(t) ~ -t 1/3 [- ln(l - p)} 2/3 (11) 

in complete agreement with the simulation results presented in figure ^ 

Before we continue, let us briefly comment on error bars and finite-size 
effects. A detailed analysis of finite-size effects in our model was performed for the 
thermodynamics in Ref. [15]. A finite-size sample contains only a finite number of 
islands of a certain size L R . As long as the number of relevant islands of size L^ p 
is large, finite-size effects on the time evolution, Fig. ^ are small and governed by 
the central-limit theorem. With increasing time, L^f increases; and the number of 
islands of this size decreases. When it becomes of order one, strong sample-to-sample 
fluctuations are expected. Finally, when L^ p becomes larger than the largest island 
on the sample, the time evolution will cross over from stretched exponential l|ll|) to 
simple exponential decay. The data in Fig. ^ show that finite-size effects do not play 
a significant role in the parameter range shown. Sample-to-sample fluctuations (as 
represented by the standard deviation of the magnetization) remain small down to 
rn = 10 -6 . Also, no deviations from stretched exponential decay are observed even at 
the smallest magnetizations shown (they do start to occur below m w 1CT 8 , though). 
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Figure 2. Time evolution of the magnetization at several temperatures from 
T = 2.8 to 3.25 (top to bottom) starting from a fully polarized state at to = to 
tmax = 2000. Parameters are as in figure The data represent averages of 25 
disorder configurations. The impurity concentrations are p = 0.2 (left panel) and 
p = 0.5 (right panel). 



We now turn to the dynamic behavior for temperatures close to but not at the 
clean critical temperature T c °. Figure [21 shows the time evolution of the magnetization 
for several temperatures from T = 2.8 to 3.25. At early times, all curves follow 
the stretched exponential found at T%. After reaching a (temperature-dependent) 
crossover time t x , the behavior qualitatively changes. For temperatures above 
the clean critical temperature, the magnetization decays faster than the stretched 
exponential. In fact, the asymptotic decay is exponential lnm ~ t/r(T), as can be 
seen from a log-linear replot of the p — 0.2 data in the left panel of figure The 
dependence of relaxation time r(T) on the temperature T is shown in the right panel 
of figure [21 Clearly, 1 /r varies linearly with T — T® in good approximation. This 
behavior can be understood as follows. For T > Tj? the correlation time of the largest 
islands does not diverge, but it is cut off by the distance from the clean critical point, 
£t ~ (T — T®)~ zv . For our infinite-range model, zv = 1. The large islands with this 
correlation time dominate the variational integral (|10|l for m(t). This leads to a simple 




Figure 3. Left: Replot of the the p = 0.2 data from figure |21 Right: Inverse 
decay time 1/t as a function of T — T° for the curves above T c °. 
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Figure 4. Double-logarithmic plot of the approach of the magnetization to the 
equilibrium value in the tail of the smeared transition for T = 2.90, 2.95 and 2.98 
(top to bottom). The impurity concentration is p = 0.5, the other parameters are 
as in figure The straight lines are fits to power laws yielding exponents of 1.78, 
2.20, and 3.31, respectively. 



exponential decay with an inverse decay time that depends linearly on T — T c °. 

The more interesting case are temperatures below the clean critical temperature 
TJ?, i.e., the tail region of the smeared transition. Here, the total magnetization 
approaches an exponentially small but non-zero value in the long-time limit, in 
agreement with Q and Ref. [15]. The approach to this stationary value is very slow. 
Fig. ^] shows m(t) — m(oo) vs. t for p — 0.5 and three different temperatures. The 
data show that the asymptotic decay is governed by a power law with a non-universal, 
temperature-dependent exponent. To understand this, we note that below T c ° the 
island correlation time has a divergence of the form !;t(T, Lr) ~ \T® — T~ AL^ 2 ^ 1 . 
If this is inserted into the variational integral , the main contribution comes from 
(finite-size) islands with L^ 2 « (T® — T)/A because they have diverging correlation 
time. Carrying out the integral leads to a power-law for m(t) — m(oo). 

The distribution of island sizes leads to a corresponding distribution of local 
correlation times Within our infinite-range model, they can be determined by 

analyzing the time evolution of the local magnetization via him x ~ —t/^ t (x). Figure 
|51 shows the resulting probability distribution P(£*) of the local correlation times at 
the clean critical point T® for p — 0.3 as determined from 10 samples of L± = 10 6 
sites. The behavior of the large-£t tail of the probability distribution follows from 
combining eq. (@J with ~ L 2 leading to InP ~ S^ 1 ^ 2 . 

4. Dynamics from general scaling arguments 

In this section, we use general scaling arguments to determine the dynamical behavior 
at the smeared transition of our disordered magnet beyond mean-field theory. The 
interesting physics in the tail of the smeared transition is local because the islands 
are effectively decoupled from each other, and the spatial correlation length remains 
finite. An appropriate quantity to study the rare region dynamics is therefore the spin 
autocorrelation function 




(12) 



X 
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Figure 5. Probability distribution of the local relaxation times at for p = 0.3, 
determined from 10 samples of size L± = 10 6 . The other parameters are as in 
figure HI The solid line is a fit of the large-£t tail to InP ~ £ t X ^ 2 . 



To determine the leading behavior of C(i) in the tail of the smeared transition we 
generalize the variational extremal statistics calculation from section [3] to the non- 
mean-field case. 

According to finite-size scaling [22], the behavior of the correlation time £ t of a 
single rare region of size Lr in the vicinity of the clean bulk critical point can be 
modelled by (for T < T c °, i.e., r < 0) 

Zt(r,L R )~L¥-W' v \r + AL- l '''\~ SP . (13) 

Here, v and z are the correlation length and dynamical exponents of a d-dimensional 
system, and v and z are the corresponding exponents of a dc-dimensional system. 

Let us first consider the time evolution of the autocorrelation function C(t) 
at the clean critical point T c °. For r = 0, the correlation time Ijl-'iH simplifies to 
£t ~ L Z R . Similar to the total magnetization in the mean-field treatment above, the 
rare region contribution to C(t) is obtained by simply summing over the exponential 
time dependencies of the individual islands with L^-dependent correlation time. Using 
Q we obtain to exponential accuracy 

C(t) ~ J dL R exp (Ljj- ln(l - p) - Dt/L R ) (14) 

where D is a constant. This integral can easily be evaluated within the saddle-point 
method. The leading long-time decay of the autocorrelation function C(t) at the clean 
critical point is given by a stretched exponential, 

lnC(t) ~-[-ln(l-p)] z /( d - L+z ) . (15) 

For T > T°, i.e. r > 0, the correlation time docs not diverge for any Lr. 
Instead, the correlation time of the large islands saturates at £t(r, Lr) ~ r~ zv for 
Lr > (r/A')~ w . The autocorrelation function C(t) can again be evaluated as an 
integral over all island contributions. For intermediate times t < t x <~ \ r \-( d ±+ z ) u ^ 
the autocorrelation function follows the stretched exponential Q15[). For times larger 
than the crossover time t x , we obtain a simple exponential decay lnC(t) ~ — t/r with 
the decay time r ~ r~ zv . Our results for T > T c ° agree with those for conventional 
Griffiths effects [7]. This is not surprising, because above Tj?, there is no difference 
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between the Griffiths and the smearing scenarios: In both cases, all rare regions are 
locally still in the disordered phase. 

This changes below the clean critical point T°. For r < 0, we repeat the saddle 
point analysis with the full expression (|13[l for the correlation length. Again, for 
intermediate times t < t x , the decay of the average density is given by the stretched 
exponential . For times larger than the crossover time t x the system realizes that 
some of the rare regions have developed static order and contribute to a non-zero 
equilibrium value of the autocorrelation function C(t). The approach of C(t) to this 
equilibrium value is dominated by finite-size islands with Lr ~ (— r/A)~ v because 
they have diverging correlation time. As a result, we obtain a power-law. 

C(t) - C(oo) ~ t~* . (16) 

The value of ip cannot be found by our methods since it depends on the neglected 
pre-exponential factors. 

The behavior of the local autocorrelation function C(x, t) = (S x (t)S x (0)) is 
spatially inhomogeneous because the local correlation time depends on the size of the 
underlying rare region. The probability distribution P(£t) of these of local (island) 
relaxation times can be obtained from and l|13fl . At the clean critical point T° we 
find 

P(6) ~ exp(-£ [- ln(l - p)} ^ /z ) (17) 
where D is a constant. 

5. Conclusions 

To summarize, we have studied the dynamic behavior of Ising magnets with planar 
defects. In these systems the magnetic phase transition is smeared because rare 
strongly coupled spatial regions independently undergo the phase transition. In the 
rare-region dominated tail of the smeared transition, the dynamics is very slow. Using 
general scaling arguments, we have shown that the spin autocorrelation function 
decays in a two stages, a stretched exponential decay at intermediate times followed by 
a power-law approach to the exponentially small stationary (equilibrium) value. We 
have illustrated these general scaling results by computer simulations of an infinite- 
range model. 

Let us briefly compare our results to the dynamics in a conventional Griffiths 
phase (as produced by uncorrelated disorder) [2,4-7]. In both cases, rare regions 
dominate the long-time dynamics. In a conventional Griffiths phase, the finite-size 
rare regions remain fluctuating for all temperatures above the dirty critical point, i.e., 
their island correlation times remain finite. As a result, the autocorrelation function 

decays like lnC(t) (lnt) d/( - d - 1 '> for Ising systems [2,4-7], and lnC(<) i 1/2 for 

Heisenberg systems [7,8]. In contrast, in the tail of a smeared transition, the effects 
of the rare regions are even stronger because individual islands can undergo the phase 
transition independently, connected with a divergent island correlation time. This 
leads to an even slower power-law decay of the spin autocorrelation function as shown 
in section 0] 

We emphasize that we have considered a purely relaxational (local) dynamics 
corresponding to model A in the Hohenberg-Halperin classification [19]. Other 
dynamic algorithms require separate investigations. For instance, model B dynamics 
globally conserves the order parameter. In this context, an interesting question is: 
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How does the interplay of the local thermodynamics of the rare regions and the global 
conservation law modify the dynamic behavior in the tail of the seared transition. 

Finally, we note that while the results here have been derived for an Ising model 
with planar defects, we expect analogous results for other disorder-smeared phase 
transitions. Indeed, in the tail of the smeared non-equilibrium phase transition of 
a contact process with extended impurities, a power-law decay of the density was 
recently found [18]. 
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